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ABSTRACT 


A mesoscale meteorological model, Gesima and the long 
range transport Acid Deposition and Oxidant model (ADOM) have 
been interfaced to obtain a nested, mesoscale air pollution 
modelling system. The system can simulate the transport, chemical 
transformations and deposition of a wide range of atmospheric 
pollutants with a horizontal resolution as fine as 5 km. In 
this study the model has been applied to simulate a high ozone 
episode in Southern Ontario, for a Toronto centered domain of 340 
x 340 km, with horizontal resolution of 20 km. The analysis of 
results focuses on areal distribution of ozone and the time 
series of ozone in selected locations, Guelph, Toronto and 
Oshawa. The results are shown to be in fairly good agreement with 


observations. 
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EXECUTIVE SUMMARY 


Two existing numerical models: a German mesoscale 
meteorological model, Gesima and the long range transport Acid 
Deposition and Oxidant Model (ADOM) have been enhanced, adapted 
and interfaced to obtain a nested modelling system for 
simulations of the transport, chemical transformations and 


deposition of atmospheric pollutants. 


The Gesima and a modified, mesoscale version of ADOM 
are ‘nested’ in the large scale ADOM domain. The role of Gesima 
is to provide the fine resolution fields of meteorological 
parameters, which for the ‘regular’, large scale/coarse 
resolution runs of ADOM are derived from the output of the 
spectral diagnostic model of the CMC at Dorval. The large scale 
ADOM provides the initial and boundary conditions for the 
mesoscale run. The interactions between Gesima and the large 
scale and mesoscale versions of ADOM are schematically 


represented in Figure l. 


To allow the nesting of Gesima in the large scale flow 
a ‘nudging’ technique has been introduced, with weak, artificial 
tendency terms forcing the mesoscale solutions towards the 


evolving large scale values. 


The role of interfacing software developed under this 
project is to extract the input data necessary for the mesoscale 
run of ADOM from the output of Gesima and the large scale ADOM. 
Software has also been developed to process other input data for 
Gesima and mesoscale ADOM, such as the topography and land use 


information. 


In order to test and validate the system, a simulation 
of a high ozone episode in Southern Ontario has been performed. 
The test run was carried out with horizontal resolution of 20 km 
in a Toronto centered area of 340 x 340 km. A period of 4 days, 


starting on August 1, 1988 was simulated. 


The analysis of results of this run focuses on the 
areal distribution of ozone and nitric oxides and their diurnal 
variations in selected locations. The horizontal cross-sections 
of the ozone and NOx concentrations, presented in Figs. 17 to 92 
show pronounced mesoscale variability, which cannot be simulated 
using the regional scale air pollution models such as the large 
scale ADOM. The modelled time series of ozone were in selected 
locations, Guelph, Toronto and Oshawa compared to the available 
observational data (see Figs. 113 - 116), and are in fairly 


good agreement with observations, especially near ground level. 
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THE DEVELOPMENT OF A LONG RANGE TRANSPORT MODEL 
WITH A NESTED FINE RESOLUTION GRID 


1. INTRODUCTION 


This report describes the development and testing of a 
nested grid, mesoscale air pollution modelling system, carried 
out under a contract with the Ontario Ministry of the 
Environment. The test run of the system concerns a high ozone 


episode over Southern Ontario in August 1988. 


The system described here is capable of simulating the 
transport, chemical transformation and deposition of a wide range 
of air pollutants, with horizontal resolution of 5 to 20 km. The 
size of the domain of simulation can vary, depending on 
resolution, from about 100 x 100 km to 500 x 500 km. The 
described test run of the system for the Toronto area, was 
carried out in a domain of 340 x 340 km, with horizontal 


resolution of 20 km. 


The resolution and the size of the domain place this system 
in the mesoscale, i.e. between the urban scale air quality models 
and the large scale models of long range transport of atmospheric 
pollutants (LRTAP). The mesoscale modelling capability is 
necessary to address the urban NOx/VOC issue with a view to 
understanding and eventually to controlling the precursor 
emissions and to quantify the subgrid scale variability of 


pollutant concentrations simulated by the LRTAP models. 


The approach was to adapt, enhance and interface two models: 
a large scale Acid Deposition and Oxidant Model (ADOM), and a 
mesoscale meteorological model ‘Gesima’. Gesima and the modified, 
mesoscale version of ADOM are then ‘nested’ in the large scale 
ADOM domain. Thus, the modelling system developed in this study 
consists of three main parts: the large scale ADOM, Gesima, and 
the mesoscale version of ADOM developed under this contract. The 
programs interfacing these three parts are also an important part 
of the system. The interactions between the system components are 


schematically represented in Figure 1. 


The original, large scale version of ADOM uses the coarse 
resolution meteorological fields from the Canadian Meteorological 
Centre (CMC) spectral diagnostic model. The mesoscale ADOM 
requires three-dimensional fields of meteorological parameters 
like horizontal and vertical wind velocity, temperature, 
humidity, turbulent diffusion parameters etc., with a fine 
horizontal resolution. Such fields are generated in this project 
by Gesima. The role of the large scale ADOM, is to provide the 


initial and boundary conditions for the mesoscale simulations. 


The development of the system and the performing of a test 
run for the Toronto area involved several tasks which are covered 


in detail in the following sections. The major tasks are: 


- modifications of the mesoscale model Gesima, allowing 
‘nested’ simulations and interfacing with mesoscale 
ADOM. 


- adaptation of ADOM, to mesoscale application. 
- development of the Gesima/mesoscale ADOM interface, to 


extract the data required by ADOM from the Gesima 
output. 


- development of the large scale ADOM/mesoscale ADOM 
interface, to extract initial and boundary conditions 
values for the mesoscale run from the output of the 
large scale ADOM. 

- assembling and processing the input data for the Gesima 
run, including the topography, land use data, and the 
meteorological data for the initial and boundary 
conditions, and for the large scale ‘nudging’ fields. 

- assembling and processing the fine resolution input 
data for the mesoscale ADOM, including the emissions 
from point and area sources and geophysical data. 

- performing the ‘production runs’ of Gesima and 


mesoscale ADOM. The large scale ADOM was not run under 
this project, the results of an earlier run were used 


instead. 

This project is a continuation of the research performed by 
M. Niewiadomski at the University of Toronto, under a 
postdoctoral fellowship from the Ontario Ministry of the 
Environment. Under that project (Niewiadomski et al., 1989, 
Niewiadomski, 1989, 1990a) Gesima was run on a smaller domain in 
the Toronto area, with a horizontal resolution of 5 km. Early 
versions of subroutines preparing the Gesima output in a form 
suitable for the mesoscale ADOM were developed, but neither the 
development of the full Gesima/ADOM and the mesoscale ADOM/large 


scale ADOM interfaces nor the running of the mesoscale ADOM were 


attempted. A separate task for this project was to rewrite all 
the JCL codes developed previously, since in the summer of 1990 
the CRAY operating system used at the Ontario Centre for Large 


Scale Computation (CLSC) changed from COS to UNICOS. 


A short description of the objectives and approach of this 
project and main features of Gesima and ADOM were presented at 
the Technology Transfer Conference in November 1990 (Niewiadomski 
1990b). The application of Gesima and its modifications developed 
for this project, especially in areas of initialization and 
‘nudging’ techniques were also presented and discussed during the 
author’s visit at GKSS in Geesthacht, Germany in April 1991. The 
‘nudging’ procedures generated some interest among the German 


scientists and may be implemented in the future at that centre. 


The following sections provide a general information about 
ADOM and Gesima models and their interfaces, describe a test run 


of the system and discuss the results. 


2. GESIMA AND ITS MODIFICATIONS 


2.1 Description 


Gesima (GEesthacht SImulation Model of the Atmosphere) is a 
sophisticated, meteorological mesoscale model developed at the 
GKSS Forschungszentrum, a research centre in Geesthacht, Germany, 


sponsored by the German federal government. 


The model solves the non-hydrostatic, anelastic equations of 
motion in terrain-following coordinates using a predictor- 
corrector scheme (for the detailed description of the dynamics 
and numerics of the model see Kapitza, 1987). The dynamic 
pressure is computed using a conjugated gradients method (Kapitza 
and Eppel, 1987). The turbulent diffusion is parameterized 
according to the mixing length approach. A sophisticated surface 
layer module is used for computing the heat and moisture balance 
at the ground level, taking into account the shortwave and 
longwave radiation, heat and moisture fluxes in soil etc. The 
model also includes a cloud module computing mixing ratios 
of cloud water, rain and cloud ice/snow (combined in one 
variable) using the ‘bulk microphysics’ approach (see Levkov et 


al, ©1969a ,.b)ic 


The version of the model, obtained from GKSS in November 
1990, had been significantly restructured and improved over the 
older version that was previously used. These changes, include 
both coding and formulational enhancements and significantly 
improved the performance and reliability of Gesima. The model has 
been recently successfully applied to study sea-breeze the in 
Northern Germany (Jacob et al, 1990) and the boundary layer 


structure over Denmark (Mengelkamp 1991). 


Information on Gesima and its modules can be found in the 
papers cited above. The full documentation of the model is 
currently being prepared and will be available from the GKSS in 
the near future. This section describes the most important 
extensions and modifications of Gesima developed under this 
project and enabling its application in the Toronto area, 


‘nesting’ in large scale fields and its interfacing with ADOM. 


2.2 Initialization 


The original Gesima is initialized by horizontally uniform 
profiles of all parameters. Since, for this project more detailed 
information from the CMC spectral model was available, those data 


were used as initial fields of wind, temperature and humidity. 


2.3 Lateral boundary conditions 


Gesima can use several types of boundary conditions, 
including fixed value at the boundary, first or second derivative 
set to zero, or the Orlanski type (where the ‘fixed’ inflow value 
is modified according to the internal flow). No ‘nesting’ 
subroutines or the time dependent boundary conditions are 


available in the original model. 


When the model is used over prolonged time period, only the 
fixed and Orlanski type conditions can accommodate evolving data 
from outside the domain, provided that the ‘fixed’ (i.e. not 
dependent on the internal conditions) inflow boundary values do 
in fact vary in time and space, according to the outer 
meteorological fields. A procedure for introducing time and 
space dependent inflow boundary values of all parameters has been 


developed under this project. 


During the initial testing of the model, it has been 
discovered that the ‘fixed’ or Orlanski type boundary conditions 
for temperature (i.e. the only types that can accommodate the 
time dependent data from outer model) cause numerical instability 
in Gesima and produce unreasonable results, especially in the 
wind field. This led to the introduction of the concept of 


‘nudging’ (see section 2.4) as an alternative way of introducing 


the time dependent information from the outer model to the 


mesoscale domain. 


Despite of the introduction of ‘nudging’ for temperature, 
humidity and wind fields, the time dependent boundary conditions 
can still be used for all of the variables except temperature. In 
the present experiment however, the time dependent boundary 


conditions were determined only for humidity and cloud 


parameters. 
2.4 Nudging 


As mentioned above, the preliminary experiments showed, that 
it is impossible, due to numerical instabilities, to introduce in 
Gesima the time dependent inflow boundary conditions for 
temperature. This led to the introduction of the Newtonian 
relaxation or "nudging" technique, not only for the temperature 


but also for wind and humidity. 


The concept of nudging (see for example Hoke and Anthes, 
1976; Stauffer and Seaman, 1990; Seaman and Cole, 1991) involves 
introduction of additional, artificial tendency terms to the 
prognostic equations of the model, which force the solution of 
the mesoscale model toward the observed or modelled large scale 
conditions. This approach allows to combine the fine scale 


variability, with the large scale tendencies which cannot be 


eae 


simulated by the mesoscale model itself, and seems to be more 
efficient and reasonable than introducing those tendencies by 


means of the time dependent boundary conditions alone. 


These artificial tendencies - the nudging terms, applied at 
each time step, are proportional to the difference between the 
values of large scale field at given grid point and the actual 
solution of the mesoscale model. The coefficient of 
proportionality - the nudging coefficient, can be a function of 


height and has values of the order of 0.0001/sec. 


3. THE MESOSCALE ADOM 


Compared to Gesima, the ADOM model is better known in Canada 
and relatively better documented (see for example Scire et al., 
1986; Venkatram et al., 1988), so this section focuses only on 
its modifications for the mesoscale application. Since the model 
is written in a general way, in principle no special, mesoscale 
version of ADOM is necessary, and its run in a mesoscale domain 
should involve only the redefining of some control parameters 
(e.g. those specifying the grid) and, of course, providing 


appropriate input data. 


However, the model has never been tested before in a nested, 
mesoscale grid, so it was expected, that many corrections or 


modifications would be required to make it behave properly in 


ige 


such a mesoscale run. Surprisingly, only a few modifications were 
necessary, including new subroutines for the time dependent 


boundary conditions. 


The ADOM input data include meteorological fields, pollutant 
emissions and their initial and boundary concentrations, geo- 
physical information, as well as some parameters controlling the 


execution of the model. 


4. SIMULATION OF THE HIGH OZONE EPISODE (August 1 = 4, 1988) 


4,1 Introduction 


This section describes the test run of the system, i.e. its 
domain, sources of the input data and the procedures of 
operation. The episode for this study, August 1 - 4, 1988 was 
chosen because of the high ozone concentrations observed at that 
time in Southern Ontario, and the availability of the results of 
the large scale ADOM simulation and the observational data. The 
model was run for a Toronto centred area. A simulation for the 


same period in the Windsor area is planned for the phase II of 


this project. 


0 


A preliminary analysis of the results of this simulation is 
provided in Section 5. It discusses some Gesima results and 
focuses on the evolution of the ozone concentrations in some 
selected locations and its areal distribution in the whole 


domain. 


4.2 The nested grid, and the domain of simulations 


The large scale ADOM is routinely run on a 33 x 33 grid 
based on a polar stereographic projection with nominal horizontal 
resolution of 127 km (true at 60° N). That domain covers most of 
eastern North America. In the Toronto area, south of 60° N, the 


actual size of this 127 km grid is about 115 km. 


In this simulation, the ‘nested’ grid of Gesima and 
mesoscale ADOM covers 9 cells of the large grid, including the 
Toronto cell and the 8 cells around it (that is the cells 15,17; 
OPEN NT, Vi peo Oe 16, Lor 17,19 215,19 to, 19 and 17; 2t9soL the 
large scale ADOM grid). Thus, the mesoscale domain covers an 
area of about 120,000 km’, with Metropolitan Toronto area located 
about 25 km south of its centre. This domain of 9 large cells is 
divided into 17 x 17 mesoscale cells, so that the horizontal 
resolution is about 20 km. A map of the area covered by the 
mesoscale domain, and its location on the large scale grid are 


shown in Figures 2 and 3, respectively. 


= qe 


The mesoscale domain applies to both Gesima and the 
mesoscale ADOM. However, while the whole domain of Gesima is 
subject of dynamic simulations (the boundary rows and columns 
numbered 0 and 18 lie outside the main 17 x 17 domain), in ADOM 
the boundary cells, with pollutant concentrations predetermined 
according to the boundary conditions, are located inside the main 
domain. Thus, the domain of simulation of ADOM is actually 
smaller than that of Gesima. This problem can be overcome in 
future simulations by redefining the domain of ADOM, (e.g. toa 
19 x 19 domain) with the boundary cells outside the domain of 


Gesima. 


The basic vertical structure of the large scale ADOM, with 
12 variable depth layers has been retained for the mesoscale 
version of the model. The same structure was adopted for the 
Gesima simulations. The tops of the ADOM layers are located at 
56.2, 135.8, 250.7, 416.3,.655.3, 1000.0, 1497.2, 2214.5 Soe 
4741.6, 6894.5, and 10000.0 metres above the ground level. The 
vertical resolution of the upper layers of ADOM is however too 
coarse for Gesima, so the 4 top layers of ADOM were split in two. 
In effect, 16 rather than 12 vertical layers are used in Gesima, 
with, tops, at. 56.2, 135.8,. 250.7, 416.3, 655.3, 1000.0, TA 
2214 .5,. 2731.8,, 3249.2,, 3995.4, 4741.6,..5818.1, 6894.5) BR 
and 10000.0 m. The Gesima/ADOM interface averages the data from 


those split layers before using them in appropriate layers of ADOM. 


Soe 


Note that the above heights of the layer tops give only the 
reference values, since both Gesima and ADOM use terrain 
following height transformations, and the real thickness and top 
height of each layer depend on the location inside the domain. 
Although the two models use different height transformations, no 
attempt has been made to include in the Gesima/ADOM interface the 
adjustment of the Gesima output to slightly different depths of 
ADOM cells. It is believed, that in relatively flat terrain of 
southern Ontario, errors resulting from this inconsistency of 
vertical grids will be negligible. The implementation of a grid 


adjustment scheme is planned for the Phase II of the project. 


4.3 Topography of the domain and the land use data 


A significant amount of time and effort was spent to develop 
software extracting the necessary data from a magnetic tape 
obtained form the Surveys and Mappings Branch of the Energy Mines 
and Resources Canada, which contained in a digital form all 
Cartographic information for a 1:250000 map. The task, while not 
impossible, proved however too cumbersome and time consuming as 
it involved extracting the coded contours of elevation from 
thousands of other cartographical features and converting the 
contours to the mean elevation in a given grid cell. The 
elevations for each of the 289 cells of the Toronto domain were 
instead subjectively read from topographical maps and manually 


entered in the data file. This procedure proved much faster than 


aie 


developing objective software. The contours of elevation in the 


mesoscale domain, obtained from these data are shown in Fig. 4. 


Land use information for both Gesima and ADOM runs was 
extracted from an EPA database. The data contain percentages of 
coverage of a given grid cell by various land use categories. 
The grid has resolution of 1/6 deg latitude and 1/4 deg 


longitude. 


4.4 Gesima simulations 


Gesima was run for a period of 96 hours starting on August 
1, 1988 at 00:00 GMT (July 31, 20:00 EDT). The time step of 45 
sec was used throughout that period. Fields of wind velocity, 
temperature and humidity were subjected to ‘nudging’ (see Section 
2.4). The nudging coefficient did not depend on height and was 
kept equal to 0.0001 s-: for all these fields except the 
temperature, for which the coefficient was increased to 


0.0005 s-1 in a part of the period of simulations. 


The fields of velocity, temperature, humidity and cloud 
water were disposed to the front end computer and briefly 
analyzed in 3 hour intervals. The ‘restart’ files were stored on 
magnetic tape in 6 hour intervals. The simulations can be 


repeated with modified input parameters starting from any of 


Ey. es 


those files. The same restart files can be used for the detailed 


analysis of any of the Gesima variables. 


With the 45 sec time step, Gesima consumed from 100 to 270 
CPU seconds of the CLSC CRAY computer, depending of the time of 
day and the meteorological situation. The total computer time of 
the 4 day simulation was about 4.5 CPU hours, i.e. roughly one 


hour per day. 


4.5 ADOM simulations 


As mentioned above, the large scale ADOM was not run under 
this project. The output files, from a previous simulation for 
the same period, were obtained from the O.M.E. These files, 
together with the output of Gesima were processed by the 
interfacing programs described in Section 3 and provided the 


input data for the mesoscale ADOM simulations. 


Since the mesoscale ADOM requires Gesima output data, the 
ADOM run started one hour later then that of Gesima, at 01:00 
GMT, on August 1, 1988 and continued for 95 hours. Thus the first 
ADOM results, after one hour of simulation, are available for 
22:00 EDT on 31 July and the last one at 20:00 EDT on August 4, 
1988. The simulations proceeded in separate runs for each day of 
the period. While the concentrations of all simulated pollutants 


stored on magnetic tape in OCON files are available for further 


ie 


analysis in hourly intervals, only 4 sets of restart files ORES 


and output files PRTOUT are available. 


The mesoscale ADOM requires much less computer time to run 
than Gesima. The simulations consumed about 1300 CPU seconds per 
day. Thus the whole 4 day GESIMA/ADOM mesoscale simulation takes 
about 6 CRAY hours plus some overhead for processing the input 


and output data. 
SE RESULTS 
5.1 Gesima results 


As mentioned above, the 3-D fields of wind velocity, 
temperature, humidity and cloud water were disposed to the front 
end computer in 3 hour intervals. The mean, maximal and minimal 
values at each level were computed from these data and compared 
to the mean, maximal and minimal values in the ‘nudging’ fields 
derived from the output of large scale CMC diagnostic model, to 
check if the modelled values stay within a reasonable range. Such 
comparison led to the mentioned above increase in the nudging 
coefficient for the temperature. The disposed 3-D fields allowed 
also to plot some two-dimensional, horizontal cross-sections, for 
wind velocity and other parameters. As an example, the wind 


velocity field at 12:00 GMT on August 1, 1988 is discussed below. 


Sie 


The contours of W-E (u) and S-N (v) components of wind 
velocity are shown in Figs. 5 and 6 for layer 1, Figs. 7 and 8 
for layer 6 and Figs. 9 and 10 for layer 10 of the model 
(elevations of the tops of each layer of the Gesima domain are 
listed in Section 4.2). For comparison, corresponding ‘nudging’ 
velocities, derived from the large scale data are shown in Figs. 
11 to 16. As can be seen, the wind field, while retaining the 
mean characteristics of the large scale winds, shows pronounced 
small scale variability, especially at the lower levels. Although 
the slightly increased wind velocity gradient over the NW shore 
of Lake Ontario can be seen in Figs. 5 and 6, the occurrence of 
lake-breezes is not clearly indicated. The 20 km horizontal 
resolution is likely too crude to simulate the local lake 


effects. 


A comprehensive analysis of the Gesima results is planned 
for the Phase II of this project, which will include numerical 
evaluations of the effects of nudging coefficients on the 


modelled fields of wind and other parameters. 


5.2 Areal distribution of ozone 


Contours of the horizontal ozone distribution in the domain, 
as modelled by mesoscale ADOM, are shown in six hour intervals in 
Figs. 17 to 64. Data from three ADOM layers are included, level 1 


denotes the first layer with mean elevation at about 25 m above 


TE 


the ground, while level 4 corresponds to about 330 m and level 7 
to about 1250 m above the ground level. Elevations of the tops of 


all ADOM layers were listed in Section 4.2. 


To explain the modelled areal distribution of ozone during 
the day and night, the concentrations of nitric oxide (NO) and 
nitrogen dioxide (NO: ) for levels 1 and 4, in 12 hour intervals 


are shown in Figs. 65 to 96. 


The concentrations of NO and NO: and some other strong 
sources peak over the area of Metropolitan Toronto and become 
transported by the wind over the area. It may be noted, that at 
the surface, the ozone is scavenged to near zero level at night 
(see Fig. 64 for 02:00 a.m. EDT, Aug. 1) in the area influenced 
by the NOx emissions, but remains over 30 ppb over the eastern 


area. 


Evolution of the meteorological situation over Southern 
Ontario and surrounding areas during the period studied is shown 
in the surface synoptic maps included at the end of this report 


(Bigs. 97 to 122). 


The strong west to east gradient of ozone concentrations 
west of Toronto, which is shown in Fig. 70 at 14:00 EDT for the 
layer 1, can be explained by the synoptic meteorological 


conditions. The leading edge of the warm air mass which contained 
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high ozone levels, moved eastward across Southern Ontario on 
Aug. 1 (see Figs. 97 to 102) reaching the Toronto area during the 
evening. Highest ozone concentrations on Aug. 1 are shown to 
occur in the south and western portion of the modelled area where 


it has been transported in the warm air mass into the Province. 


5.3 Diurnal variations 


The diurnal variations of the modelled ozone concentration 
is easier to analyze by concentration vs time graphics. Such 
plots have been prepared for a few selected locations and are 


discussed below. 


To show the ozone level as determined by the large scale 
ADOM, an output has been produced for a cell on the western 
boundary of the modelled area, where the pollutant concentrations 
follow the boundary conditions, determined according to the large 


scale ADOM results. 


It may be noted (see Fig. 113), that the large scale ADOM 
has produced the first peak of surface ozone on Aug. 1 at 18:00 
EDT (78 ppb), the second peak at 15:00 EDT on Aug. 2 (72 ppb), a 
minor peak (50 ppb) 04:00 Aug. 3, which accompanied a major peak 
(128 ppb) aloft, another peak (100 ppb) on Aug.3 at 20:00 EDT and 


a peak of 70 ppb on Aug. 4 at 15:00 EDT. The timing of some of 
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the peak ozone levels at the surface at the western edge of the 


area are difficult to explain from meteorology during the period. 


For a location near Guelph, the model output at three 
elevations for ozone, nitric oxide and nitrogen dioxide is shown 
in Figs. 114, 115 and 116. It may be noted that the first peak in 
surface ozone occurs at 21:00 EDT on Aug. 1, possibly with the 
arrival of the warm polluted air, ozone then peaks at the more 
expected times of 15:00 EDT on Aug. 2, 14:00 EDT on Aug. 3 and 
maintains persisting high concentrations, near 80 ppb, from 13:00 
to the end of simulations at 20:00 EDT on Aug. 4, 1988. It may be 
noted, that ozone at elevated levels peaked at almost the same 
times as at the surface, except during the early morning of 
Aug. 4 when ozone peaked aloft at the time when ozone at the 


surface was at a minimum. 


5.4 Comparison with observations 


In order to validate the ozone levels as computed by the 
ADOM/Gesima approach, the discussed above concentrations at the 
lowest model cell in the Guelph area are compared in Fig. 117 to 
the levels measured at the AQI station in Guelph. Similar 
comparisons are shown in Fig. 118 for the surface of the Toronto 
cell of the model and measurements at 26 Breadalbane in Toronto, 
in Fig. 119 for the level 4 elevation in Toronto and the 


measurements at the CN Tower, and in Fig. 120 for the surface 
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modelled concentrations in the Oshawa area and the observations 
at the AQI station in Oshawa. Note, that due to the 
unavailability of the observed data for some hours, the periods 
shown in Figs. 117 to 120 may differ from those discussed in 


Section 5.3 and between each other. 


Analyzing these figures, one should keep in mind, that while 
the measurements were taken at specified points, the modelled 
concentrations represent averages from a mesoscale ADOM cell of 
the size 20 x 20 km and the depth of 56 m for the ‘level’ 1 
and 165 m for ‘level’ 4. Even so, the modelled and observed ozone 
concentrations are in fairly good agreement, especially for the 


surface locations in Toronto and Oshawa (Figs. 118, 120). 


The times of maximum and minimum modelled concentrations 
coincide with those observed at the surface in Toronto and Oshawa 
(except for the secondary peaks in observed concentrations during 
the early morning of Aug 4, both in Toronto and Oshawa, which are 
not reflected in the model results). For the Guelph location 
(Fig. 117) the modelled peak occurs 3 hours after the observed 
one on Aug. 1, coincides on Aug. 2 and roughly on Aug. 4 and 


advances the observed peak by about 3 hours on Aug. 3. 


The values of the modelled ozone concentrations are also in 
fairly good agreement with the observed ones, although the model 


does not predict the highest measured concentrations (e.g. 112 


wo te 


ppb observed in Toronto on Aug. 2 at 14:00 EDT versus 68 ppb 
modelled). These discrepancies can be partly attributed to a 
significant area and depth represented by an ADOM cell, compared 


to the point measurements. 


The upper level data, both modelled and observed, have a 
much less regular diurnal variability (see Fig. 119). The 
modelled concentrations roughly agree with the observed ones, the 
biggest discrepancies being the peaks in the model results at 
03:00 EDT on Aug. 2 and at 08:00 EDT on Aug 4 during the 


secondary and primary minima in observations. 


6. CONCLUSIONS 


A mesoscale modelling system for transport, chemical 
transformations and deposition of atmospheric pollutant has been 


designed, programmed and tested in one case study over Southern 


Ontario. 


The system is capable to simulate the distribution of 
various air pollutants with the spatial resolution of 20 km or 
less. It can detect local effects in this scale, that are 


impossible to simulate with large scale models like the original 


version of ADOM. 


399 = 


The preliminary analysis of the simulation of the early 
August 1988 high ozone episode over Southern Ontario, showed 


reasonably good agreement of the model results with observations. 
More comprehensive analysis of results of several case 
studies, and further improvements of the system are planned for 
the Phase II of this project. 
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Figure 1: Components and scheme of operation of the 
nested mesoscale modelling system. 
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Figure 2: The mesoscale domain. 





Figure 3: Domain of the large scale ADOM and location 
of the mesoscale domain. 
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Figure 4: The mesoscale domain, contours of elevation. 
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